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Quantum interference (QI) effects in molecular junctions may be used to obtain large 
thermoelectric responses. We study the electrical conductance G and the thermoelec¬ 
tric response of a series of molecules featuring a quinoid core using density functional 
theory (DFT), as well as a semi-empirical interacting model Hamiltonian describing 
the TT-system of the molecule which we treat in the GW approximation. Molecules 
with a quinoid type structure are shown to have two distinct destructive QI features 
close to the frontier orbital energies. These manifest themselves as two dips in the 
transmission, that remain separated, even when either electron donating or withdraw¬ 
ing side groups are added. We hnd that the position of the dips in the transmission 
and the frontier molecular levels can be chemically controlled by varying the electron 
donating or withdrawing character of the side groups as well as the conjugation length 
inside the molecule. This feature results in a very high thermoelectric power factor 
S'^G and figure of merit ZT, where S is the Seebeck coefficient, making quinoid type 
molecules potential candidates for efficient thermoelectric devices. 


I. INTRODUCTION 

Calculations suggest that quantum inter¬ 
ference (QI) effects in molecules give rise to 
a large tunability of the electron transport 
properties of molecular junctions. These ef¬ 
fects induce a strong variation of the trans¬ 
mission with energy, which is favourable 
for thermoelectricity, where an electric cur¬ 
rent or potential difference develops in re- 
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sponse to a temperature difference across 
the moleculei^-. Well studied QI molecu¬ 
lar units usually involve either cross conju¬ 
gation and/or meta coupled phenyl units^^-. 
Molecules such as anthraquinone, have re¬ 
cently been shown to exhibit destructive QI 
effects in a molecular junctiort^'i^. A sim¬ 
ple nearest- neighbour tight-binding model of 
the TT system has been used to rationalise the 
QI effects. Moreover, schemes have been de¬ 
veloped to make predictions based on sim- 
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pie graphical The sign of the ther¬ 

mopower has been suggested to provide infor¬ 
mation of whether the transport is mainly via 
occupied -or unoccupied molecular states^. 
To obtain a strong thermoelectric response, 
the destructive QI feature needs to be close 
to the Fermi level. A sensible handle to con¬ 
trol the position of QI features is offered by 
tuning the chemistry of the binding groups or 
changing the electronegativity of substituent 
side groups, which has been shown theoret¬ 
ically as well as experimentally to influence 
the thermopower by changing the molecular 
frontier levels relative to the Fermi level of 
the metal electrodes^^— . 

Cruciform ohgo(phenylene ethynylene)s 
(OPEs) with a conjugated and extended 
tetrathiafulvalene (TTF) donor moiety have 
recently been synthesised^'^ and shown to 
be redox active as well as having interest¬ 
ing spin properties in the Coulomb block¬ 
ade regime^^. Cruciform type molecules have 
also been synthesised with substituent side 
groups, such as TTF, dithiofulvalene (DTF) 
and atomic oxygen, forming cross-conjugated 
OPEs. Such structures can be referred to as 
quinoid, since the central core corresponds 
to a quinone, with the substituents replac¬ 
ing the oxygen atoms, see Figure [I^ and 
Figure |5l They may exist as zwitterions 
with mixed substituents, an electron donor 
on one side and an electron accepting group 


on the other side. This suggests the possibil¬ 
ity of electric field induced switching between 
the conjugated (high conductance) and cross- 
conjugated (low conductance) state. While 
previous theoretical studies have found that 
the thermoelectric response of molecules may 
be greatly enhanced by QI effects, they in¬ 
volved either rather long molecular wires, 
radicals or metal complexes without anchor- 
mg group3^““. 

Here we explore the electron transport 
properties of cruciform molecules with a 
quinoid type structure by varying the elec¬ 
tron donating (ED) and electron withdrawing 
(EW) character of substituent side groups. 
We study the low bias conductance G as 
well as the ability of the molecules to con¬ 
vert thermal energy into electric energy by 
applying a temperature difference across the 
molecular junction. We hnd, based on den¬ 
sity functional theory (DFT) calculations, 
that quinoid structures show two charac¬ 
teristic destructive QI (DQI) features, one 
near the highest occupied molecular orbital 
(HOMO) level and one near the lowest unoc¬ 
cupied molecular orbital (LUMO) level. The 
DQI features result not only in a high ther¬ 
mopower (Seebeck coefficient) S', but more 
importantly from the point of view of techno¬ 
logical application, a high power factor (PF), 
S‘^G. The PF is related to the electrical work 
that can be extracted from a thermoelectric 
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device and is the quantity that determines 
the dimensionless hgure of merit, ZT = 
S'^GT/(Kph + Kei), when the phonon ther¬ 
mal conductance Kph dominates over the elec¬ 
tron thermal conductance We hnd 

that a semi-empirical model for the vr sys¬ 
tem, treated in the GW approximatioii^^^— , 
gives similar maximal PFs as the DFT cal¬ 
culations. We note that a simple nearest- 
neighbour tight-binding or Hiickel model 
does not capture the split interference fea¬ 
ture. Quinoid type structures typically yield 
a maximal PF an order of magnitude higher 
than similar molecules showing DQI in the 
centre of the HOMO-LUMO gap, such as 
meta coupled benzene or simple acyclic cross 
conjugated molecules. 

II. METHOD 

We use DFT as implemented in the 
GPAW and ADF codes to provide an quan¬ 
tum chemical description of the charge 
transport through the molecular junction 
system22i2S. We also use a more ap¬ 
proximate density functional tight-binding 
(DFTB) method^. The three DFT based 
methods allow us to get an estimate for 
how sensitive the results are to the par¬ 
ticular implementation details. Molecules 
were optimised in the gas phase using the 
LDA exchange correlation (xc) functional in 


ADF and the PBE^^ xc functional in GPAW. 
Structures used in the DFTB method were 
relaxed using the B3LYP xc functional. In 
all calculations the molecules were attached 
to the FGG hollow site of Au(lll) with a Au- 
S bond length of 2.5 A (1.83 A above the sur¬ 
face). In GPAW, the scattering region super¬ 
cell was modelled modelled using 3-4 atomic 
An layers on both side of the molecule. The 
number of surface layer atoms varies between 
4x4 and 6x6 depending on the size the 
molecule and periodic boundary conditions 
where used in the transverse directions. The 
2D Brillouin zone was samples using 4x4 
fc-points. In ADF the extended molecular re¬ 
gion includes 3x3 An atoms in the surface 
and no periodic boundary conditions are used 
in the transverse directions. The An atoms 
were frozen in the bulk lattice structure us¬ 
ing the DFT derived lattice constant (PBE: 
a = 4.18 A, LDA: a = 4.08 A). 

For calculations based on GW, a semi- 
empirical model Hamiltonian based on the 
Pariser-Parr-Pople (PPP) type for describ¬ 
ing the TT-system is used^. We use the 
Ohno parametrisatioii^ with long range two- 
electron interactions with the parameter U = 
10.0 eV describing the Goulomb repulsion 
and nearest neighbour hopping element of 
t = —2.5 eV. On site energies are taken rela¬ 
tive to carbon, which we set to ec=0. For the 
GW calculations we use a wide band approx- 
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imation for the leads. More details about the III. RESULTS 


GW method can be found in Reference 
28. 
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Transport properties are for all 
methods calculated using the Lan- 
dauer Buttiker transmission formula 
expressed in terms of Green’s functions 
T{e) = Tr[G'"(£)r^(£)G'“(£)r^(£)], where 

r“ = — E“) is given in terms of the 

lead a self-energy Eq. We calculate the 
isothermal conductance in the zero bias 
voltage limit as G = Go / {—n'p{e,T))T{e)de, 
where n'p{e,T) is the derivative of the Fermi 
function with respect to energy. T is the 
temperature and Go = 2e^//i is the unit of 
quantum conductance, where h and e are 
Planck’s constant and the electronic charge, 
respectively. We calculate the thermopower 
in the limit of low temperature drop across 
the molecule from the transmission as 
^ = {2e/hT) J e{-n'p{e,T))T{e)de/G. The 
power factor (PF) is calculated as S^G and 
we use room temperature T = 300 K, see 
Appendix |B] for more details. We note that 
we do not use the Sommerfeld expansion 
expression for the calculations of the ther¬ 
mopower, since the transmissions functions 
we consider may have structure on the scale 
of keT^. 


We begin by considering a simple model 
for a quinoid structure, which captures 
the essential physical mechanism responsi¬ 
ble for interesting thermoelectric properties 
of quinoids. The schematic structure of a 
transport junction with a quinoid molecule 
sandwiched between two leads i shown in Fig¬ 
ure [T^. The connection of the central unit to 
the leads are indicated by the dashed lines. 
Ri and i ?2 denote substituent side groups. 
In Figure [Tb we show the PF (S'^G) calcu¬ 
lated with GW for a PPP model Hamiltonian 
where Ri and R 2 are taken as GH 2 groups, 
as indicated by the inset. The model in¬ 
cludes a single orbital per carbon atom 
all with equal on site energies. For compar¬ 
ison we also show the PF obtained for ben¬ 
zene coupled in para and meta position. The 
maximal PF, within the HOMO-LUMO gap, 
for the quinoid structure is seen to reach a 
value that is an order of magnitude higher 
than both meta and para coupled benzene. 
In fact the PF for the quinoid is compa¬ 
rable to the maximal PF that can be ob¬ 
tained from a Lorentzian transmission line 
shape, ~ 0.9 k|/h., by using a width of about 
1.1 keT, see Figure fTTh in Appendix iBl 

The PFs can be analysed in terms of the 
transmission and thermopower as a function 
of energy, shown in Figure [It and [T]i, re- 
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FIG. 1. (a) Sketch of a quinoid molecule sand¬ 
wiched between metal leads. Ri and R 2 denote 
the substituent side groups, (b) Powerfactor 
S'^G for a quinoid structure and benzene con¬ 
nected in ortho and meta position. The insets 
shows the chemical structures of the molecules. 
The results are obtained using a PPP model with 
GW. (c) and (d) show the transmission and ther¬ 
mopower as a function of energy. 

spectively. The transmissions for both the 
quinoid and the meta coupled benzene ring 
is seen to be highly suppressed, as compared 
with para coupled benzene, for an extended 
energy region around the Fermi level. This 
is a result of DQI, as we will discuss in more 
detail below. DQI also introduces a strong 
variation of the transmission with energy 
which in turn yields a high thermopower. 
This can be seen by considering the linear 
response Sommerfeld expansion expression 


for the thermopower S{e) oc clg log[T(e)p^. 
The meta-coupled benzene ring has a single 
dip in the transmission at the centre of the 
HOMO-LUMO gap, while the quinoid struc¬ 
ture shows two transmission dips - one dip is 
close to the HOMO energy while the other 
dip is close to the LUMO energy. While 
the thermopower for meta-coupled benzene 
reaches high values close to the Fermi level, 
the PF ends up being low because the trans¬ 
mission, and thus G, is very low here. The 
quinoid structure on the other hand, has the 
DQI feature near a transmission resonance 
giving both a high thermopower and trans¬ 
mission within the same energy region. This 
combination is responsible for the high PF for 
the quinoid structure. The para-coupled ben¬ 
zene shows a relatively low thermopower, but 
relatively high transmission and thus ends up 
having a PF very similar to the meta-coupled 
benzene. Finally we note that the quinoid 
structure has a smaller HOMO-LUMO gap, 
which is related to the longer conjugation 
length as compared with benzene, he. it is 
a conhnement effect. The maximal PF value 
is obtained at an energy located between the 
HOMO resonance and the dip, which is about 
0.2 eV away form the HOMO resonance for 
all three molecules. 

We now analyse how the thermoelectric 
properties can be improved by adding side 
groups to the quinoid backbone. The aim 
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is to move the DQI feature close to a trans¬ 
mission resonance (HOMO or LUMO) which 
leads to a peak in the PF, and to shift that 
peak towards the metal’s Fermi energy. 

To probe how the two DQI features de¬ 
pend on the choice of the side groups, we 
show in Figure [2^ and [2b how varying the on¬ 
site energies of the side groups in the model 
Hamiltonian affects the transmission and PF, 
respectively. The side group site energies are 
0.0 eV for both Ri and R 2 in the top panel, - 
4.0 eV for both Ri and R 2 in the middle panel 
and -4.0 and 4.0 eV for Ri and R 2 in the 
lower panel. The negative/positive site ener¬ 
gies are meant to mimic the effect of EW/ED 
side groups. We note that ED groups (EDG) 
and EW groups (EWG) will in general also 
tend to shift the on-site energies of the ring 
through inductive and resonance effects by 
perturbing the electron density within the 
ring in the a and tt- system, respectively. The 
perturbation of the electron density through 
the cr-system is caused by polarisation and is 
not included in the PPP model. The change 
in electron density will in turn affect the elec¬ 
trostatic potential and can be seen as a way 
of chemically gating the molecular backbone. 

The position of the transmission dips is 
seen to shift with the side group energies 
but the two split QI features are quite ro¬ 
bust and stay within the HOMO-LUMO gap. 
It can be seen that a transmission dip can 




FIG. 2. (a) Transmissions for the quinoid type 

structure with on-site energies (Fii, E 2 ) of the 
{Rii R 2 ) side groups of (0.0, 0.0) top, (-4.0,-4.0) 
middle and (-4.0,-1-4.0) bottom in units of eV. 
The vertical dashed lines indicates the HOMO 
and LUMO levels, (b) Same as in (a) but for 
the power factor, S‘^G. 

be moved closer to the HOMO resonance 
by lowering the side group energies, which 
leads to an increase in the corresponding PF. 
This can be rationalised from the approxi¬ 
mate expression S (x log[T(e)] for the ther¬ 
mopower, that is, the slope of the transmis¬ 
sion as function of energy on a logarithmic 
scale. The slope of the transmission is seen 
to increase in Figure |2b when on-site ener¬ 
gies are varied such that the transmission 
anti-resonance moves closer to the molecu¬ 
lar frontier energy levels, i.e. the transmis¬ 
sion has to change from a value of 1 at reso¬ 
nance to 0 at the antiresonance over a smaller 
energy region. The calculated maximal PF 
within the HOMO-LUMO gap are high in 
all cases, and reach a value of ~ 1.1 ke^/h 
(=316 fW/K^). A similar value was obtained 
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by optimising a general two level model in 
Reference y. We therefore expect that the 
EDGs and EWGs can not only be used to 
tune the position of the interference dips, 
but also to enhance the PE. To experimen¬ 
tally obtain high PFs, molecules are needed 
where either the HOMO or the LUMO level 
are close to the Fermi level. This can be 
chemically controlled by the binding (anchor¬ 
ing) groups^^— , or by mechanical tuning^, 
as well as by the EW/ED side groups. Before 
exploring a more detailed quantum chemical 
description of the molecules shown in Fig¬ 
ure [SI we analyse the two split interference 
features in more detail for the model of the 
quinoid core structure. 

The appearance of the QI features can 
be analysed in a number of ways ranging 
from non-spanning nodes to interfering path¬ 
ways through Feynman path integrals and 
the phases of molecular orbitals^'^^— . Here 
we rationalise the QI interference features in 
terms of molecular orbitals and correspond¬ 
ing molecular energies. This has the advan¬ 
tage that in the weak molecule lead coupling 
limit this analysis can be generalised for the 
correlated electron case by considering so- 
called Dyson orbitals and energies^. 

We argue that the central quinoid core 
structure is responsible for the two transmis¬ 
sion nodes. This can be seen either based 
on graphical rules developed to predict QI 


in Hiickel models!^ or by systematically re¬ 
ducing the Kohn-Sham Hamiltonian atom 
by atom as will be discussed below. We 
note that while a nearest neighbour Hiickel 
model predicts degenerate transmission anti¬ 
resonances for equal site energiesl^, both the 
GW and the DFT calculations give two well 
separated anti-resonances for quinoid struc¬ 
tures. A similar effect was observed and for¬ 
mulated in terms of long range hopping ele¬ 
ments for ac;w;lic cross-conjugated molecules 
in Reference 


For simplicity we proceed with an effec¬ 
tive non-interacting model where we use the 
Hiickel molecular orbitals but quasi-particle 
energies taken from the GW calculations. 
The use of quasi-particle energies instead 
of Hiickel molecular yields split interference 
features^. The two central DQI features 
of the quinoid structure, see Figure [I^,c, 
can be accounted for by considering the 
four molecular orbitals closest to the Fermi 
level (HOMO-1, HOMO, LUMO, LUMO+1). 
This may be seen by hrst considering the 
condition for the transmission going to zero, 
T{e) (X |G/r(e)P = 0. Here the Green’s func¬ 
tion Gir describes the probability amplitude 
for an electron or hole propagating through 
the molecule to a site I from a site r on the 
molecule, which are connected to the left and 
right lead, respectively. Then, by express¬ 
ing Gir in terms of the molecular (uncoupled) 
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Green’s function as Gir{e) = / D{e)^, 

we see that a necessary condition for a trans¬ 
mission zero is that is zero for some 

energy. The denominator in the expression 
for Gir is given by D = (1 - " 

In the spectral representation 
may be expressed in terms of molecular 
orbitals {V'n} and energies {£„} as 


Gr‘(E) = E 

n 


e - En + ir]' 


( 1 ) 


where 77 is a positive inhnitesimah The ele¬ 
ments {l\'4>n) and {r\'il>n) give the amplitude 
of the n’th molecular orbital on site I and r, 
respectively. 

From Eq. ([1]) we infer that the relative 
sign of the HOMO- and LUMO amplitudes 
on the sites / and r determine the parity 
of the number of transmission zeros within 
the HOMO-LUMO gap. This is illustrated 
in Figure |3^ and [3 ]d, where we for a few 
representative cases sketch the real part for 
of within the HOMO-LUMO gap with 
an even and odd number of zero crossings, 
respectively. The contribution to Re[G™°^] 
from the HOMO and LUMO diverges when 
approached from within the HOMO-LUMO 
gap towards either -|-cxo or — 00 , with the sign 
determined by the residue (=(/|'0n)('0n|r))). 
All other orbitals only contribute with a h- 
nite value to Re[G™°'(£)] within the HOMO- 
LUMO gap. This means that the function 
Re[G™°\e)] can in general be drawn as a 


continuous line connecting the divergence at 
^HOMO and ^lumoj with the detailed shape 
determined by all HOMOs and LUMOs. If 
the residues of the HOMO and LUMO have 
different signs we are in the situation corre¬ 
sponding to Figure [3^, where we see that any 
continuous line is forced to cross zero and 
even number (0,2,4,...) of times. With four 
orbitals, as in our model, only 0 or 2 zero 
crossings are possible. 

If the the HOMO and LUMO residues 
have the same sign, we are in the situation 
shown in Figure [3b, where any continuous 
line is forced to cross zero an odd (1, 3, 5,...) 
number of times. From this a quick assess¬ 
ment and classihcation of interference fea¬ 
tures can be made in terms of the HOMO 
and LUMO orbitals only: Constructive inter¬ 
ference between HOMO and LUMO results in 
an even number of transmission zeros while 
for destructive interference the number of ze¬ 
ros is odd. 

To simplify the analysis even further we 
now turn to the Coulson-Rushbrooke (CR) 
pairing theorem^ for alternating hydrocar¬ 
bons, i.e. molecules which can be viewed 
as bipartite. The pairing theorem relies 
on a particle-hole symmetry and states that 
the molecular energies come in pairs as 
^HOMO-n = ^LUMO+n for u = 0,1, 2,... . The 
molecular orbitals of such a pair are identical 
except for a sign change on one of the sub- 





FIG. 3. (a) Sketch of a few representative cases 
of Re[G“°*] from HOMO and LUMO orbitals 
having different relative sign on sites connected 
to leads, (b) same as (a) but for HOMO and 
LUMO orbitals with same sign on sites con¬ 
nected to leads. 

lattices, i.e. every other atom, which means 
that the sign of the residues of the orbitals in 
a CR pair can be predicted without any cal¬ 
culations. In short, if leads are connected to 
the same sub-lattice of a molecule the orbitals 
in a CR pair, they interfere destructively, giv¬ 
ing an odd number transmission zeros. For 
leads connected at sites belonging to different 
sub-lattices, the two orbitals in a CR pair in¬ 
terfere constructively, giving an even number 
of transmission zeros. In FigureIDD, we sketch 
the four important quinoid frontier molecular 
orbitals, with the weight on a site represented 
by the size of the circle and the sign by the 
color. The dashed lines indicate the sites / 
and r, for which leads are connected. The 
starred sites constitute one sub lattice of the 
molecule. We see that the paired orbitals, 
namely the HOMO and LUMO and also the 


HOMO-1 and LUMO-l-1, indeed follow the 
behaviour predicted by the CR theorem. Due 
to the close relation between orbitals in a CR 
pair, we hnd it advantageous to hrst sum up 
the contribution from a CR pair in Eq. ([T]) 
and subsequently sum up all pairs. 

In order to enable destructive interference 
in the case where CR-paired orbitals inter¬ 
fere constructively, two or more sets of such 
paired orbitals are needed, i.e. a minimum of 
four molecular orbitals. We show the paired 
HOMO-LUMO (blue line) and the paired 
HOMO-1 and LUMO-l-1 (red line) contribu¬ 
tions to Re[G™°^] for the quinoid in Figure 
0^. Here the vertical dashed lines indicate 
the positions of the quasi-particle levels. It 
can be seen that the HOMO- LUMO pair 
adds up constructively due to the different 
signs arising from the product of the molec¬ 
ular orbital weights on the I and r sites - 
a consequence of the leads being connected 
to different sub- lattices. In the same way 
the HOMO-1 - LUMO-l-1 pair interferes con¬ 
structively. The destructive interference ob¬ 
served in the total transmission arise when 
we consider the sum of the pair contribu¬ 
tions (green line), which is seen to cross zero 
twice in accordance with the general classih- 
cation based on the amplitudes of the HOMO 
and LUMO orbitals on the sites connected to 
leads. The transmission nodes are a conse¬ 
quence of the HOMO and LUMO pair hav- 
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FIG. 4. (a) Contributions to the real part 

of the molecular Green’s function for the 

quinoid structure; HOMO and LUMO (blue 
curve) pair and the HOMO-1 and LUMO-l-1 pair 
(red curve). The vertical dashed lines indicate 
the positions of the poles, arising from the molec¬ 
ular energy levels, (b) The molecular orbitals 
are represented with the weight on the different 
sites by the size of the circles and the sign by 
the colour, (c) and (d), same as (a) but para 
and meta connected benzene, respectively, (e) 
same as (b) but for benzene para. 

ing a smaller weight on the I and r sites 
than the HOMO-1 and LUMO-l-1 pair. If we 
imagine increasing the HOMO and LUMO 
orbital weight on the lead sites, then at some 
point the transmission nodes will disappear, 
althongh the transmission will still be snp- 


pressed. The DQI is for the qninoid therefore 
not only a result of the phases of the molec¬ 
ular orbitals. Similarly the relative position 
of the HOMO- and HOMO-1 level (and the 
LUMO- and LUMO-l-1 level) also influence 
the appearance of transmission zeros. 

It is interesting to compare the quinoid to 
benzene in para position, where the orbitals 
making up a OR pair interfere constructively 
just as for the quinoid, see Figure St andH^. 
However, in contrast to the quinoid, adding 
up the two pairs now leads to constructive 
interference. This is a result of a different 
phase relations between the two OR pairs 
for benzene. For meta coupled benzene, the 
CR pairing theorem ensures that the two or¬ 
bitals making up a CR pair always interfere 
destructively, resulting in a zero crossing in 
Re[G“°\e)] at the Fermi level, see Figure IDi. 

As already mentioned, a nearest neigh¬ 
bour Hiickel or tight-binding model, predicts 
a doubly degenerate node, and does not cap¬ 
ture the characteristic split interference fea¬ 
ture of the quinoid structurel^*^. In terms of 
the analysis presented here, which was based 
on Eq. the doubly degenerate transmis¬ 
sion node can be interpreted graphically as 
in Figure UK- Re[G™°^] (blue line) may be 
shifted up along the y-axis by a change in 
the molecular energy levels (or a change in 
orbital weight on / and r). At some point, the 
two zero crossings coincide at the Fermi level. 
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This is exactly what happens if one nses the 
tight-binding molecular orbitals and energies 
to construct The splitting can be ob¬ 

tained by adding a scissors operator to the 
H iickel Hamiltonian to change the HOMO- 
LUMO gap. The scissors operator will in 
general have long range hopping matrix el¬ 
ements, which has been shown to split inter¬ 
ference features in some cases^. 

From the four frontier molecular orbitals 
in Figure |Dd, we can estimate the effect of 
changing the on- site energies on the side 
groups. The HOMO-1 and LUMO-l-1 have 
(almost) zero weight on the side group sites 
and a perturbation of the on-site energies, 
Vsg, will therefore not change the HOMO-1 
and LUMO-l-1 orbitals and energies to low¬ 
est order since Vsg\a) 0 for a G {HOMO- 
1, LUMO+1}. Thus, the HOMO- and 
LUMO level shift a lot, but the HOMO-1 and 
LUMO-l-1 do not. This can be used to give 
a graphical interpretation in Figure [2{i. A 
change of Ae of the side group on-site ener¬ 
gies will result in a shift of the HOMO- and 
LUMO levels of Ae, while the HOMO-1- and 
LUMO-l-1 levels stay fixed. This will result 
in a shift of the position of the two dips in 
the same direction as the HOMO and LUMO 
level. 

To explore the robustness as well as the 
tunability of the DQI features with respect to 
the EW/ED character of the substituent side 



FIG. 5. Molecules M1-M5 with a quinoid type 
structure consisting of a central cross-conjugated 
phenyl ring. Ml contains CH 2 side groups, M2 
contains a TTF electron donating unit, M3 con¬ 
tains a TTF unit with esters, M4 is a possible 
zwitter-ion with an electron donating DTF unit 
with ester groups on one side, and an oxygen 
atom on the other side. Finally M5 contains two 
electron withdrawing oxygen atoms. 

groups we have performed quantum chemical 
calculations based on DFT for the molecules 
M1-M5 shown in Figure [5l 

Before discussing molecules M1-M5 in 
more detail, we give a comparison between 
the transmission calculated with the different 
DFT methods in Figure |6t for Ml. The two 
split DQI features discussed for the simple 
quinoid model above are clearly visible and 
we emphasise the very good agreement be¬ 
tween the GPAW and ADF calculations de¬ 
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spite small variations in geometry and the use 
of different xc functionals. We note that the 
agreement between the methods, especially 
the position of the Fermi level, is less good 
for M2-M5. However, the side group induced 
changes and trends are not affected and we 
therefore proceed discussing the results ob¬ 
tained with GPAW. 

We stress that we have used flat Au(lll) 
surfaces for simplicity and not investigated 
the effect of using different lead structures, 
such as a small pyramids or add atoms on a 
flat surface. While the use of different lead 
structures may certainly change the conduc¬ 
tance and the thermoelectric properties we 
do not think it will change the trends and 
conclusions as these are based on a property 
of the molecule itself. 

To make a more direct connection to the 
minimal Pz model for the quinoid core struc¬ 
ture, we have performed calculations start¬ 
ing from the molecule Ml and then system¬ 
atically removed atoms from its arms. The 
transmissions in Figure Eb are obtained by 
only considering the Pz subspace of the Kohn- 
Sham Hamiltonian of the full calculation for 
successive truncations of the arms of the 
molecule Ml, while keeping the central cross- 
conjugated quinoid unit unchanged. The 
different transmissions correspond to the 
molecules denoted Mla-Mld shown in Fig¬ 
ure Et. Wide band leads are attached to the 


atoms connected with a dashed line. The 
position of the transmission nodes are seen 
to be quite insensitive to the arms of the 
molecule attached to the central quinoid core 
structure. Only when leads are attached di¬ 
rectly to the quinoid core structure do we see 
a non-negligible change. We stress that in 
all cases the qualitative shape of the trans¬ 
mission function remains intact and in good 
agreement with the GW results in Figure [2b. 

In Figure [7b and [Tb we show the trans¬ 
mission and PF calculated using GPAW (see 
Appendix [A] for the corresponding ADF and 
DFTB transmission results) for the molecules 
M1-M5. The vertical lines indicate the posi¬ 
tion of the molecular frontier energy levels 
obtained by diagonalising the molecular sub¬ 
space (including sulfur). 

For all molecules, except M5, we see clear 
signs of the two split DQI features. Inter¬ 
estingly, M5 only shows a single transmission 
anti-resonance, which according to the dis¬ 
cussion above based on Eq. ([T]) implies that 
the HOMO and LUMO orbital in this case 
interfere destructively. Taking molecule Ml 
as the baseline, we observe that the HOMO 
energy level is closer to the Fermi level than 
the LUMO energy level. This is typical for 
molecules bonded to An via thiols^. The 
similarity of the shape of the transmission 
with the results for the simple quinoid model 
(see Figure [2] a) is striking. We see two DQI 
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FIG. 6. (a) Transmission as a function of en¬ 
ergy for molecule Ml calculated with GPAW, 
ADF and DFTB. (b) Transmissions within the 
Pz subspace but for truncated versions of Ml as 
denoted in (c). Wide band leads are attached at 
the Pz orbitals connected to ” An” with a dashed 
line in (c). (d) Frontier molecular orbitals for 

Ml. 


features, one located near the HOMO and 
the other near the LUMO level. The PF 
for Ml is shown in the upper most panel 
in Figure [7 ]d and has a maximal value of 
~ 1.5 within the HOMO-LUMO gap. 

This is larger than, but comparable to our 
GW result for the quinoid model and about 
60% higher than the maximal value of a sin¬ 
gle level model with a Lorentzian line shape. 
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FIG. 7. (a) Transmissions for molecule M1-M5 
calculated with DFT (GPAW). (b) Power factor 
for M1-M5. The vertical colored lines indicate 
the position of the frontier energy levels of the 
molecule (HOMO: red, LUMO: blue). 


see Appendix [Bl Interestingly, the maximal 
PF for Ml is even larger than the opti¬ 
mized value obtained for a general two-level 
non-interacting model^. The high PF ob¬ 
tained 0.75 eV below the Fermi level is a re¬ 
sult of the transmission anti-resonance being 
close to a resonance such that both the ther¬ 
mopower and conductance are high at the 
same energy, as already discussed for the sim¬ 
ple quinoid model. In Figure |6]i we show 
the frontier molecular orbitals for Ml which 
are responsible for the two split transmission 
nodes. The orbitals can be compared with 
the simple model in Figure |Dd and are seen to 
have the same essential features: The HOMO 
and LUMO orbitals have smaller weight at 
the connecting sites than the HOMO-1 and 
LUMO-l-1 orbitals, and the orbitals in a CR 
pair interfere constructively. Also the sym- 
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metry of the orbitals with respect to the 
transport direction are the same. The essen¬ 
tial mechanism behind the two DQI featnres 
in the DFT calculations is therefore captured 
by the simple PPP semi-empirical vr-system 
model for the quinoid core structure. The PF 
for Ml at the Fermi level is rather low, sug¬ 
gesting that the band alignment of the molec¬ 
ular levels with the Fermi level of gold is not 
ideal in terms of possible thermoelectric tech¬ 
nological applications. 

Molecule M2 contains two electron- 
rich five membered rings containing sulfur 
(DTFs) that have an ED character. These 
units are therefore expected to push the 
molecular levels up in energy as compared 
with Ml without the ED rings. Indeed, we 
see that the HOMO -and LUMO levels are 
moved up in energy with the HOMO level 
ending up very close to the Fermi level. We 
note that with the GPAW code, the HOMO 
level ends up just below the Fermi level, 
while it ends up above the Fermi level with 
the ADF code. We have verified that this 
is not a basis set issue, by comparing to 
real space grid calculations performed with 
GPAW. Based on this we speculate that this 
difference is due to use of a cluster model 
instead of a periodic structure for the leads 
in ADF. The HOMO-LUMO gap of M2 is 
considerably reduced as compared with Ml. 
This is an effect of the increased conjuga¬ 


tion length of the TTF side groups in M2 
compared with the GH 2 side groups in Ml. 
However, the two split DQI features are still 
intact and within the HOMO-LUMO gap. 
While the maximal PF is similar to Ml, the 
PF at the Fermi level is now relatively high, 
with a value of 0.4 k|/e (~ 100 fW/K^). 

Molecule M3 is similar to M2 but with es¬ 
ters added to the DTF substituents which 
may have a benehcial effect on the solubil¬ 
ity of the molecule. The addition of esters 
alters the transmission spectra slightly near 
the Fermi level as compared to M2. How¬ 
ever, there are now unoccupied transmission 
resonances at higher energies, such that the 
high-energy anti- resonance now falls out¬ 
side the HOMO-LUMO gap. These trans¬ 
mission resonances do not reach a value of 1 
and originate from two orbitals, LUMO and 
LUMO-l-1, asymmetrically coupled to leads 
and with most of the weight localised on the 
DTF-|-ester side group. The PF at the Fermi 
level is again rather high and comparable to 
what we found for the M2 molecule. 

The M4 molecule is an interesting 
molecule and a possible zwitterion. It con¬ 
tains an EDG (DTF) as one side group and 
an EWG (atomic oxygen) as the second side 
group. This gives the molecule a large dipole 
moment. The transmission as a function of 
energy shows the two DQI transmission fea¬ 
tures. The HOMO is coupled rather asym- 
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FIG. 8. (a) vr-system contribution to the trans¬ 
mission for M5. The vertical dashed lines indi¬ 
cates the positions of the molecular energy levels, 
(b) Frontier molecular orbitals for M5. 

metrically which results in a transmission res¬ 
onance with a value below 1. We see that 
the opposite character of the two side groups 
tends to cancel the overall shift of the levels, 
in agreement with the PPP model studied 
above, see Figure [2^. 

Molecule M5 has a real quinone core. A 
similar molecule, namely the anthraquinone 
was recently shown to exhibit DQl features 
in the dl/dV-i^. The two electronegative 
oxygen atoms are EW which is reflected by 
the molecular levels being drawn down in en¬ 
ergy. The LUMO is now closer to Ep than 
the HOMO. There is clearly only one trans¬ 
mission dip within the HOMO-LUMO gap. 
In Figure |8^ we show the tt system contribu¬ 
tion to the transmission for M5, from which 


it is clear that a second interference feature 
is present at —2.5 eV below the Fermi level. 
However, since this is not within the HOMO- 
LUMO gap it is not relevant. By inspecting 
the frontier molecular orbitals in Figure |Hb 
we see that that while the LUMO is the same 
as for Ml, the HOMO is not. The symme¬ 
try of the HOMO and LUMO orbitals now 
implies destructive interference within the 
HOMO-LUMO gap. From the analysis based 
on Eq. m, we expect an odd number of 
transmission zeros within the HOMO-LUMO 
gap in this case. This is consistent with the 
observation of one transmission node. The 
transmission node is positioned close to the 
Fermi level which induces a strong variation 
of the transmission function and thus a large 
thermopower. However, the transmission is 
suppressed, and thus the resulting PF at the 
Fermi level is rather low. 

To investigate the general thermoelectric 
performance of the molecules we show in Fig¬ 
ure inti- and inb the dimensionless hgure of 
merit ZT as a function of energy for molecule 
Ml and M2, respectively. To probe the de¬ 
pendence on the phonon thermal conduc¬ 
tance, Kph, we show the results for four dif¬ 
ferent values ranging from 0 to 100 pW/K. 
Thermal conductance values in the range 10- 
100 pW/K have recently been calculated for 
similar OPES molecules^^. Measurements on 
n-Alkanes suggest a similar range of values. 
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TABLE I. ZT and PF for the molecules MI- 




FIG. 9. Figure of merit for Ml and M2 is shown 
in (a) and (b), respectively. The different lines 
are obtained by using different values for the 
phonon thermal conductance Uph, as indicated 
in the legends. The inset shows the maximal 
value of ZT within the HOMO-LUMO gap as 
function of the phonon thermal conductance. 


M5 obtained using DFT-PBE. ZT is for Kpj^ = 
0,10 pW/K. PF is in units of (~ 

288fW/K^). 


Molecule (ZT)max 

(ZT)ep 

(PF)max 

Ml 

4.1, 

2.7 

0.014, 0.001 

1.5 

M2 

4.5, 

2.8 

4.2, 1.8 

1.2 

M3 

4.1, 

2.7 

3.5, 1.3 

1.1 

M4 

5.5, 

2.5 

0.032, 0.004 

0.7 

M5 

1.8, 

1.5 

0.706, 0.005 

1.0 


We see that maximal ZT, denoted (ZT)max, is 
almost identical for the two molecules. This 
is especially clear when comparing the in¬ 
sets which show (ZT)max as a function of 
Kph- For comparison the maximal ZT which 
can be obtained using a thermal conductance 
of Kph = 10 pW/K for a single level model 
with a Lorentzian line shape is about 2.4, see 
Appendix |Bl However, in this case a rather 
narrow resonance is needed with a with of 
about O.SkeT. The dependence of ZT on en¬ 
ergy, suggests that molecular levels need to 
be aligned with the Fermi level with a pre¬ 
cision of about ~ 0.1 eV. For molecule M2 
the alignment predicted by DFT (GPAW) is 
relatively good giving high values for ZT at 
the Fermi level. We have collected (ZT)max 
and ZT evaluated at the Fermi level, (ZT)ep, 
and (PF)max in Table HI Here we clearly see 


that all the quinoids show promising thermo¬ 
electric properties when the maximal values 
are considered. In terms of the level align¬ 
ment predicted by DFT molecules M2 and 
M3 look promising both in terms of power 
production and efficiency. We note that the 
level prediction by DFT may not be very ac¬ 
curate, however, we believe that the trend in 
the level position induced by the substituent 
side group is a robust feature. Also, we have 
only considered one type of binding structure 
between the sulfur anchoring group and the 
gold surface. Both the level positions and the 
broadening of levels may be sensitive to the 
details of the gold-sulfur interface and this 
will in turn affect the calculated thermoelec¬ 
tric properties. However, if it is mainly the 
broadening of the levels that are affected by 
the binding geometry, then we do not expect 
large changes in the calculated properties. 
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This is because the destructive QI forces the 
transmission to change from 1 to 0 between a 
frontier orbital resonance and a nearby anti¬ 
resonance irrespectively of the broadening. 

IV. CONCLUSIONS 

We have shown, based on quantum chemi¬ 
cal calculations using DFT as well as GW for 
an interacting semi-empirical Hamiltonian, 
that molecules with a quinoid topology may 
show very high power factors and ZT values, 
which suggests a high power generation per 
molecule and good efficiency, respectively. 

The good thermoelectric properties were 
found to originate from a particular DQI fea¬ 
ture of quinoid type structures, namely two 
split interference features within the HOMO- 
LUMO gap, resulting in a transmission anti¬ 
resonance lying close to a molecular reso¬ 
nance. This feature was shown to involve 
the four nearest frontier molecular orbitals. 
The split interference feature was found to 
be rather robust and the position of the reso¬ 
nances and anti-resonances was shown to be 
highly tuneable by the ED/EW nature of side 
group substituents. For the five molecules 
studied, only the real quinone core breaks 
the two split interference feature, only hav¬ 
ing a single transmission node in the HOMO- 
LUMO gap. 

DFT calculations showed that the max¬ 


imal power factor as well as ZT obtain¬ 
able within the HOMO-LUMO gap was only 
weakly dependent on the chemical nature 
of the side group. However, by varying 
the ED/EW character of the side group the 
power factor and ZT values evaluated at 
the Fermi level can be tuned. Three differ¬ 
ent DFT based methods predicted the same 
trends, however, the exact position of the mid 
HOMO- LUMO gap relative to the Fermi 
level could differ by up to ~ 0.3 eV. The 
maximal ZT values were predicted by DFT 
to vary from 1 to 3, for phonon thermal con¬ 
ductances in the range 10 to 100 pW/K. 


The high power factors we predict may 
be affected by inelastic transport channels, 
such as those arising from the interaction of 
electrons with phonons. Here the current 
from inelastic processes may be signihcant, 
since the HOMO resonance is close to the 
anti-resonance. An incoming electron with 
an energy at the anti-resonance may emit a 
phonon with an energy that brings it close to 
the HOMO resonance and thereby bypass the 
anti- resonance. We have neglected such pro¬ 
cesses, but they appear to be import aspects 
for further studies of thermoelectric proper¬ 
ties of molecules relying on anti-resonances 
close to resonances. 
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FIG. 10. Transmission calculated using GPAW 
(black line), ADF (red line), DFTB (green line). 
(a)-(e) corresponds to molecule M1-M5. 

Appendix B: Thermoelectrics 


The transport coefficients relevant for 
thermoelectricity is written as 


Appendix A: Transmission comparison 

We compare in Figure [TOl the transmission 
calculated for molecule M1-M5 using GPAW, 
ADF and DFTB. We observe that the GPAW 
and ADF results are in overall good agree¬ 
ment, and that all three methods predicts 
the same qualitative behaviour and trends 
in accordance with the electron withdraw¬ 
ing and donating character of the substituent 
side groups. 


G = e^Lo 

(Bl) 

eT Lq 

(B2) 

ILl 

Ke = Ko - - — 

1 Lq 

(B3) 

= Ko- TGS^ 

(B4) 

^0 = 

(B5) 

in terms of the function 



where np{e) = (exp((e — fi)/kBT)) + 1)“^ is 
the Fermi-Dirac distribution function. The 
thermal electronic conductance ko is for zero 
chemical potential drop and related to the 
thermal conductance at zero electric current 
by Ke = K .0 — TGS"^. 
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The dimensionless figure of merit which 
may be used to characterise the performance 
of a thermoelectric device is given by 

TGS^ 

ZT = 


Ke T f^ph 

^ / 1 - Ke/tio \ 

^ph \ 1 ^e/ ^ph ) 

< —, (B8) 

f^ph 


though this results in a ZT below the MS 
bound. 

Single level model. We now consider a 
single level coupled to wide band leads for 
which the transmission takes the a simple 
Lorentzian form 




(B9) 


which we have rewritten in terms of the ther¬ 
mal conductances in the second line. Since 
all thermal conductances are larger or equal 
to zero, we see immediately that the expres¬ 
sion in the parentheses is always smaller or 
equal to 1. This means that Ko/i^ph is an 
upper bound to ZT for a given transmission 
function and phonon thermal conductance. 
This upper bound, reached when Ke = 0, 
is sometimes referred to as the Mahan-Sofo 
(MS) bound^. MS showed that the only 
a Dirac delta function as the transmission 
function gives ZT = Ko/^e, i-e. the upper 
bound. However, when concerned with sin¬ 
gle molecule junctions, then not only is the 
integral of the transmission bounded but also 
the value at any energy. The transmission is 
usually smaller than ~ 1 at the relevant en¬ 
ergies within the HOMO-LUMO gap. This 
has the important consequence, that in the 
case of a infinitesimally narrow transmission 
resonance kq —^ 0. As we shall see below, 
typically a finite width of the order of keT 
results in the highest figure of merit, even 


where T gives the broadening and Sa is the 
level position. We assume a temperature of 
T = 300K unless otherwise stated. In Figure 
[nk we show the maximal power factor as a 
function of T. The level position relative to 
the Fermi level yielding the maximal power 
factor for a given F is shown in Figure fTTb . 
We have used a temperature of lOOK (blue 
line) and 300K (dashed red line) which gives 
the same results. Note that energy is in units 
of ksT. The largest value for the maximal 
power factor is ~ 0.9 ke^/h = 258/hF/JF^ 
obtained for F ps 1.1 keT and the level po¬ 
sitioned about 2.9 keT away fro the Fermi 
level. 

We show in Figure ITTh the maximal value 
of ZT, denoted as (ZT)max, as a function F 
for different phonon thermal conductances, 
Kph, as indicated in the legends. The dashed 
lines shows the MS upper bounds, and we 
see that only in the limit of small F does the 
calculated (ZT)max approach the MS bounds. 
However, in this limit ZT goes to zero. The 
largest value of (ZT)max is obtained for a res- 
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FIG. 11. (a) Maximal power factor as a function 
of the broadening F for T = lOOif (blue line) and 
T = SOOiC (red dashed line), (b) Level position 
relative to the Fermi energy giving the maximal 
power factor in (a) for a given F. (c) Maximal 
figure of merit as a function of F for T = 300 K. 
The labels indicate the values of the contribution 
to the thermal conductance from phonons. The 
dashed lines show the MS upper bound, (d) The 
level position giving the maximal figure of merit 
in (c) for a given F. 

onances with a with of about kBT/2 and the 
level energy positioned about 2.5 k^T away 
from the Fermi level, see Figure fTTH. We ob- 
serbe that while (ZT)max diverges in the limit 
of zero broadening for Kph = 0 pW/K the 
highest possible value for a finite Kph, quickly 
decreases. 
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